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We assess the capabilities of magnetic Resonant Inelastic X-ray Scattering (RIXS) at the O K 
edge in undoped cuprates by taking La2Cu04 as a benchmark case, based on a series of RIXS 
measurements that we present here. By combining the experimental results with basic theory we 
point out the fingerprints of bimagnon in the O K edge RIXS spectra. These are a dominant peak 
around 450 meV, the almost complete absence of dispersion both with tt and a polarization and 
the almost constant intensity vs. the transferred momentum with a polarization. This behavior is 
quite different from Cu L3 edge RIXS giving a strongly dispersing bimagnon tending to zero at the 
center of the Brillouin zone. This is clearly shown by RIXS measurements at the Cu L3 edge that 
we present. The Cu L3 bimagnon spectra and those at Cu K edge - both from the literature and 
from our data - however, have the same shape. These similarities and differences are understood 
in terms of different sampling of the bimagnon continuum. This panorama points out the unique 
possibilities offered by O K RIXS in the study of magnetic excitations in cuprates near the center 
of the BZ. 

PACS numbers: 78.70. Ck, 78.70.En, 74.72.Cj, 75.30.Ds 



I. INTRODUCTION 



In the debate on high T c superconductivity the role of 
the magnetic excitations/fluctuations and the interplay 
between charge and spin degrees of freedom continue to 
be among the central issues. 1 ' 2 Quite recently this debate 
has been boosted by the application of Resonant Inelastic 
X-ray Scattering (RIXS) 3,4 which extended the experi- 
mental possibilities offered by well established methods 
as angle resolved photoemission 5 and neutrons. 6 On the 
experimental side a breakthrough has been the demon- 
stration that it is possible to study magnetic excitations 
(namely bimagnon) in cuprates with RIXS. 7 This new 
approach cross fertilizes especially with neutron spec- 
troscopy actively used in the study of magnons (see recent 
work in Refs. 8-11) and expands the experimental pos- 
sibilities; in fact generally speaking x-rays allow a wide 
exploration of g-space with easy detection of high en- 
ergy excitations and require very tiny amount of sample 
material as opposite to neutrons. On the other hand the 
major drawbacks are the energy resolution which is much 
worse than with neutrons and the competition with flu- 
orescence which becomes overwhelming in good metals 
making the application of RIXS to magnetic excitations 
very difficult or even impossible. For these reasons RIXS 
is particularly suitable in (quasi) bidimensional cuprates 
since the low energy scale is expanded by the high su- 
perexchange and the RIXS signal is found up to high 



dopings covering, to say the least, a great fraction of the 
dopings corresponding to the superconducting dome in 
the phase diagram. 

Successively it has been shown experimentally 12 and 
explained theoretically 13,14 that RIXS at the Cu L3 
edge gives direct access to single magnon excitations. 
Ref. 12 shows moreover that in the underdoped regime 
magnon-like excitations reminiscent of the excitations 
in the parent compound are seen up to high energies 
(typically 350 meV). The opportunities offered by RIXS 
have been recently exploited with Cu L3 RIXS in the 
study of the YBa2Cu306+^ family showing that damped 
magnons (i.e. paramagnons) are present in all cases stud- 
ied therein spanning from underdoped to slightly over- 
doped regime. 15 Moreover, Ref. 15 shows that these para- 
magnons are compatible with very high T c superconduc- 
tivity supported by spin fluctuations. As far as higher 
order spin excitations are concerned a mixture of single 
magnon and bimagnon has been observed with Cu L3 
RIXS, 16 while it has been demonstrated that RIXS at the 
Cu K edge gives experimental access to bimagnon 7,17 in 
agreement with the theoretical work of Refs. 18-21. Thus 
with a suitable choice of the method it is possible to study 
either single magnon and high order odd excitations or 
bimagnon and high order even excitations offering alto- 
gether a detailed panorama of spin excitations in bidi- 
mensional cuprates. Needless to say this possibility is of 
paramount importance in the debate on superconductiv- 
ity. 
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The topic of the present paper and of the successive 
one is mainly RIXS at the oxygen K edge respectively 
in a benchmark parent compound, the La2Cu04, and in 
the derived high T c systems, the La2- x Sr x Cu04 fam- 
ily. In the panorama sketched above it might seem that 
RIXS at the O K edge has nothing or little to add to the 
above scenario. Here we demonstrate that this is not the 
case and that O K RIXS offers unique opportunities. It 
is in fact possible to map the even order spin excitations 
dominated by an almost non-dispersing bimagnon with 
a peak around 450 meV in the low q range (up to about 
40% of the BZ boundary). We show that this region of 
the parameter space is not accessible to Cu L 3 and Cu K 
RIXS while it is known to be, to say the least, unfavorable 
to neutron scattering. On top of that traditional Raman 
spectroscopy is limited to the case of zero momentum 
transfer. The present work is timely because of the rapid 
evolution of the field and needed because very little on 
spin excitations seen with RIXS at the oxygen K edge can 
be found in the literature where a systematic treatment 
is lacking: in this connection we mention the bimagnon 
feature measured with O K RIXS in S^CuC^C^, re- 
ported for the first time in 2002 by Harada et al? 2 and 
more recently by Guarise et al. 23 

Our results on O K RIXS are presented in two twin 
papers. Here (paper I) we assess the method and by com- 
paring with RIXS at the other edges we obtain a com- 
prehensive overview of bimagnon spectroscopy. To this 
end we present also original measurements at the Cu K 
and L 3 edges. In the spirit of a mostly phenomenological 
paper we take advantage of a basic theoretical approach 
based on Ref. 19. A byproduct is the way of disentan- 
gling mono and bimagnons in Cu L3 spectra. In paper II 
we exploit the peculiarities of O K RIXS to study a typi- 
cal superconductor family i.e. the doped La2- x Sr x Cu04. 
Thus the present paper identifies the fingerprints of even 
order spin excitation seen at oxygen K edge. 

The present paper is organized as follows. We give 
the experimental details for RIXS at the three edges in 
Sec. II. The O K data are given in Sec. Ill, where they 
are presented (subsection III A) and discusssed (subsec- 
tion IIIB). Next, in Sec. IV data measured at Cu L3 edge 
and Cu K edge are given with a comparison between 
O K and Cu L3 (in subsection IV A) and between Cu 
K and Cu L3 (in subsection IV B). The conclusions are 
drawn in Sec. V. The paper is supplemented by two ap- 
pendices: Appendix A presents the method to extract 
the bimagnon excitation at the Cu L3 edge; Appendix B 
deals with multiphonons contributions at the O K edge. 



II. EXPERIMENTAL 

The La2Cu04 (LCO) sample used in soft x-ray RIXS 
at the O K and Cu L 3 edges was a 100 nm film grown by 
pulsed laser deposition on (001) SrTi03, and is the same 
used in Ref. 12. The sample measured with hard x-rays 
at the Cu K edge was a single crystal grown by travelling 




FIG. 1. (Color online), a) The geometry of the RIXS exper- 
iment: the sample is in green with its crystallographic axis 
in blue, while the scattering plane is in light blue. The red 
and the black arrows indicate the incoming photon polariza- 
tions, respectively parallel (a) and perpendicular (ty) to the 
(a,b) plane, b) The cuprate Brillouin zone. The dashed line 
marks the magnetic BZ boundary. The colored area highlight 
the q region of the 1 st BZ accessible at the O K edge in the 
present mounting (green), the Cu L3 edge (lilac) and the Cu 
K edge (orange). In c) and d) the two different experimental 
geometries correspond to the maximum momentum transfer 
condition (5 — —50°) and the zero parallel momentum trans- 
fer condition (5 = 0°). 



floating zone method. 24 In both cases the surface was the 
basal (a, b) plane with the c axis in the scattering plane 
(see Fig. 1). 

The O K RIXS measurements (photon energy ~ 
530 eV) were done at the beamline ID08 of the ESRF 
equipped with the AXES spectrometer; 25 the measure- 
ments were taken in part before and in part after an 
intervention on the equipment. In the first phase the 
AXES spectrometer was looking (with 110° scattering 
angle) to the sample before the entrance slit of the 
Dragon monochromator of the beamline and the incident 
monochromatic light was given by the dedicated Polifemo 
monochromator. 26 In the second phase AXES was sitting 
at the end of the beamline past the Dragon monochroma- 
tor with a sagittal refocusing done by a toroidal mirror. 
In this last case the scattering angle is 130°. The com- 
bined energy resolution was adjusted to the same value 
(around 150 meV) in both setups and the new installa- 
tion resulted in a 3-4 times higher count rate. In the 
present paper the spectra taken with the second set up 
are used unless otherwise stated. The average time per 
spectrum with the Dragon set-up was typically 4 hours 
for a polarization i. e. with the electric vector in the (a, b) 
plane. 
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The Cu L 3 RIXS measurements were done at the 
ADRESS beamline 27 of the Swiss Light Source equipped 
with the SAXES spectrometer, 28 with a scattering angle 
of 130°, and a combined energy resolution of 130 meV 
at 930 eV. Each spectrum took 30 min. A summary 
of the geometrical configuration used for experiments in 
the soft x-rays is given in Fig. 1(a). The total momentum 
transfer is Q = ki n — k ou t- Being the sample basically 
bidimensional, the interesting component of the momen- 
tum is given by q, the projection of Q onto the basal 
(a, b) plane. This is controlled by choosing the scattering 
angle a and by rotating the sample with 0. Hereafter the 
angular rotation is indicated by the angle S between the 
sample surface normal and the bisector of the scattering 
angle. Thus q oc sinS (see Fig. l(c-d)). 

The measurements at Cu K edge were done at ID 16 
of the ESRF equipped with the RIXS spectrometer of 
Refs. 29 and 30, in Rowland circle geometry. The ex- 
perimental configuration was chosen according to Ref. 7: 
we used only tt type polarization in order to maximize 
the projection of the incoming photon wave vector onto 
the c axis (at this edge the bimagnon is observed only 
if the incoming polarization has a component along that 
axis 7 ) and a scattering angle close to 90° to reduce the 
elastic contribution. The spectra were measured near the 
second Brillouin zone (BZ) indexed (2,2,11) in order to 
avoid diffraction peaks and to have the reduced momen- 
tum transfer q completely defined in the (a, b) plane. The 
combined energy resolution was adjusted to 170 meV for 
this experiment. The average acquisition time per spec- 
trum was 8 hours. In the following, the Cu K spectra will 
be presented versus the reduced momentum q reported 
to the first BZ. All measurements at all edges were taken 
at 20 K. 



III. THE BIMAGNON AT O K EDGE 

The RIXS study presented here is carried out with ex- 
citation in the threshold region of the oxygen K edge 
absorption spectrum. The absorption at this edge has 
been studied extensively, 31 33 because it sheds light into 
the valence states hybridized between O 2p and Cu 3d 
states. This is the origin of the pre-edge peak shown by 
the arrow (labeled Eul) in Fig. 2a and is characteristic of 
undoped cuprates. This excitation is used in the present 
RIXS work, which is indeed a resonant study since, as we 
will show, there are spectral features suppressed or dras- 
tically reduced with the excitation energy above the peak. 
The general shape of the emission spectrum is given in 
Fig. 2b. As it is known most of the intensity belongs to 
the resonant fluorescence 34 ' 35 which covers the range be- 
tween about 3 eV and 10 eV energy loss when the sample 
is excited as in Fig. 2. Moreover the broad peak around 
2 eV is the superposition of ^-excitations 36 on the Cu 
site and of charge-transfer exciton of the kind of a Zhang- 
Rice singlet. 22,37 At even smaller energy loss i.e. in the 
Mid Infrared (MI) range (blue box in Fig. 2b expanded 




energy loss (eV) energy loss (eV) 



FIG. 2. (Color online), a) O K absorption, b) Emission 
spectrum of La2Cu04 measured with a polarization at S = 
0°. The excitation energy, Eul, is marked by the arrow in 
a), c) The low energy expansion of the spectrum shows the 
bimagnon feature. The red arrow points to its peak position, 
around 450 meV. 



in Fig. 2c), the spectrum shows the feature we are inter- 
ested in. In this region we see a dominant peak around 
450 meV whose dependence on excitation energy, inci- 
dent polarization and transferred momentum is presented 
hereafter. We will show that all these results are consis- 
tent with the assignment of this peak to a bimagnon ex- 
citation together with higher order spin excitations of the 
same symmetry i.e. with an even number of spin flips. 
When no misunderstanding is possible we will simplify 
the language by calling conventionally bimagnon the ex- 
perimental result including also the higher orders. 



A. Experimental results in the Mid Infrared 

The main experimental results are summarized in 
Fig. 3. The panel 3a shows the polarization dependence 
of the absorption; note that the spectra are consistent 
with the literature. 31 33 The excitation energy depen- 
dence is given in panel 3b showing the resonant behavior 
of the MI feature at the peak excitation Eul with a dras- 
tic decrease at excitation energy Eu2, above the peak 
(data taken at 110° scattering angle). The excitation 
energies are given by the vertical bars in panel 3a. 

The polarization and the momentum dependence of 
RIXS are shown in panels 3c-f where the black spectra 
give the tt incident polarization and the red spectra the a 
one (data taken at 130° scattering angle). These are raw 
data normalized to the acquisition time. These results 
show two characteristic aspects. First the dispersion with 
the transferred moment q is, if any, very small both along 
the direction from T to (1,0) and from T to (1,1); this is 
summarized in panels (g) and (h) including also the re- 
sults at 110° scattering angle (q=(0.3,0) along (1,0)). A 
weak but visible dispersion is seen only with tt polariza- 
tion along the (1,0) direction. Second the spectral weight 
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FIG. 3. (Color online), a) XAS spectrum of La2Cu04 for a (red) and tt (black) polarizations (5 — —45°) measured in 
fluorescence yield mode, b) On resonance (red open dots) and out of resonance (red filled diamond) RIXS spectra measured 
respectively at Eul and Eu2 excitation energies, marked in a). The blue line is a smoothing of the inelastic part, as a guide 
to the eye. c-f) O K RIXS spectra measured at Eul as a function of momentum and polarization (scattering angle 130°). c) 
and e) RIXS spectra along the (1,0) direction respectively for a (red) and tt polarization (black); d) and f) refer to the (1,1) 
direction. Spectra for q ^ (filled symbols) are vertically offset with respect to the spectra for g=0 (open symbols). All the 
spectra have been normalized to the acquisition time, g-h) Summary of the bimagnon dispersion respectively along the (1,0) 
and along the (1,1) direction. The data are expressed as a function of q x , the component of the in-plane momentum q along 
the (1,0) direction. Data points for g=(0.3,0) and (0.21,0.21) refer to the spectra measured at 110° scattering angle. 



of the MI feature is basically determined by the absorp- 
tion, i.e. the excitation channel with tt polarization has a 
typical loss of intensity at grazing incidence (high q). In 
this case the electric vector has a small component in the 
(a,b) plane so that the probability of exciting the system 
decreases, being the ground state of x 2 — y 2 type with 
the holes in the (a,b) plane. 



B. Discussion 

The assignment of the MI excitations was the object of 
long standing discussions on the relative weights of po- 
laronic effects and of magnetic excitations; 38,39 now the 
situation is settled at least in RIXS at Cu L3 and K edges 
with a compelling evidence of the assignment to magnetic 
excitations. As a matter of fact at Cu K edge the low- 
est magnetic excitation is the bimagnon since there is no 
spin-orbit so that only an even number of spin flips be- 
sides no spin flip is allowed; in this connection the small 
spin orbit in the valence state of 3d systems has negligi- 
ble effect. On the other hand at Cu L3 edge the core hole 
spin-orbit is large and the spin is not a constant of the 



motion; thus the single magnon is clearly seen with its 
characteristic dispersion in agreement with the neutron 
results. Obviously this framework and the analogy with 
Cu K edge indirect excitation of bimagnon is a strong 
support to the assignment of the MI feature in O K RIXS 
to a bimagnon. We stress that in both cases there is no 
spin-orbit. The observation in Sr2Cu02Cl2 at O K edge 
and assignment to bimagnon presented at a single q value 
by Guarise et al 23 is essentially based on these qualita- 
tive arguments. Well before this work a faint feature 
has been seen and correctly assigned in Sr2Cu02Cl2 by 
Harada et al 22 but the experimental basis was at that 
time very limited due to technical constraints. 

A first strong support to this assignment is the compar- 
ison of the RIXS results with the density of states (DOS) 
of the bimagnon; to this end the non-interacting results 
of Fig. 4 (reproduced from Ref. 19) are sufficient. The O 
K RIXS process covers a limited region of q around the T 
point. In the figure the two green boxes give the regions 
accessible at the O K edge due to the limited available 
momentum. In this range the DOS has a prominent in- 
tensity at high energy and a weak intensity at very low 
energy with a continuum in between. Unless dramatic 
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FIG. 4. (Color online). The calculated two-magnon density of 
states for nearest-neighbor interaction and without magnon- 
magnon interaction, as reported in Ref. 19. The green rectan- 
gles highlight the BZ portions investigated at the O K edge. 
The energy scale on the left is in units of J, while the one on 
the right is in eV (J=135 meV 41 ). 



effect of the matrix elements comes into play, the sam- 
pling of the bimagnon would bring to a dominant fea- 
ture at high energy. This would be somewhat lower than 
4J, where J is the superexchange. Being this calcula- 
tion without interaction, it is expected that the energy 
scale is overestimated with respect to the experiment, as 
it will be discussed in the following. With this sampling 
the absence of dispersion or a small dispersion are a quite 
natural consequence. 

Going into more details the assignment to the bi- 
magnon is supported also by a basic theoretical model 
already used in the indirect bimagnon excitation at Cu 
K edge. 19 This is based on the linear term in a series 
of powers of J/T representing the scattering operator (T 
is the core level lifetime broadening). With respect to 
Ref. 40 one introduces the oxygen ions and the hopping 
from O 2p to Cu 3d. Also in this case, in order to un- 
derstand the type of sampling, the inclusion of magnon- 
magnon interaction is not needed. The theory confirms 
that one is sampling the upper intense ridge of the DOS. 
This is shown by the comprehensive view of the theo- 
retical cross sections corrected for self-absorption and in- 
strumental resolution given in Fig. 5 by color maps in nice 
qualitative agreement with the experimental results. In 
fact basically no dispersion is seen with a polarization 
while a tiny dispersion might be compatible with the tt 
excitation along (1,0). Also the trend of the intensities is 
well reproduced by the theory. As a matter of fact along 
the (1,0) direction with a polarization also the small in- 
crease of the intensity seen in the experiment is obtained. 
This is due to the reduction of self-absorption at grazing 
incidence and is an effect of about 20%; note that this 
is the only case in which the self-absorption correction 
really matters. 
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FIG. 5. (Color online). RIXS theoretical cross sections calcu- 
lated as in Ref. 40 with 130° scattering angle. The experimen- 
tal broadening and self-absorption correction are included. 
The columns refer to the inspected BZ directions, (1,0) and 
(1,1), while the rows refer to the polarizations a and tt. The 
red vertical bars indicate the instrumental resolution, ~ J. 
As in Fig. 3, q x is the component of the in-plane momentum 
q along the (1,0) direction. 



Having assessed the O K edge RIXS it is interesting 
to compare with the other edges at Cu. This is done in 
the next section. 



IV. COMPARISON WITH THE OTHER EDGES 

In order to obtain a general view of bimagnon spec- 
troscopy, another piece of experimental information is 
missing at present. This is the bimagnon information 
one obtains from Cu L3 edge RIXS. Indeed the work at 
Cu L3 edge is emerging as a powerful tool in the study of 
single magnon. As a consequence, up to now minor at- 
tention has been given to the information on bimagnon 
basically left dormant in the raw data. Here we present 
results on the bimagnon at the Cu L3 edge in a wider 
context i.e. the comparison between the spectra of bi- 
magnon at three edges i.e. O Cu L3 and Cu K. This 
is an interesting problem since at each q value, the bi- 
magnon density of states is represented by a continuum 
so that different sampling are typical of RIXS at each 
edge and there is no a priori reason why they should be 
equal. To make this program viable we need a way to ex- 
tract the bimagnon information from the L3 data where 
it is superimposed to single magnon. More precisely we 
need a method to separate the contributions from odd 
and even number of spin flips. Also this problem is ad- 
dressed in the present paper and for better readability 
this methodological contribution is presented separately 
in Appendix A. 
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FIG. 6. (Color online), a) The bimagnon extracted from the 
Cu L3 RIXS spectra along the (1,0) direction with a polariza- 
tion, b) Cu L3 XAS of LCO. c) The dispersion as seen with O 
K edge (red squares) and Cu L3 (black squares) for the (1,0) 
direction. The gray arrows point out the correspondence to 
the theoretical two-magnon DOS along (1,0), in d). 19 



A. Comparison between O K and Cu L3 RIXS 

This comparison is done with q along the direction 
from T to (1,0) and the main results are summarized 
in Fig. 6. In panel 6a we give the bimagnon spectra 
extracted from the L3 data with the method given in 
App. A; the excitation is at the L3 peak in the absorption 
as specified in panel 6b. Note that the bimagnon spectra 
extracted as in Appendix A have the same spectral func- 
tion independently of the polarization which determines 
simply a scaling factor. The results of panel 6c show a 
clear positive dispersion at L3 as opposite to the small 
negative dispersion at higher energies seen in O K data 
and already discussed. This striking difference is under- 
stood once again in terms of different sampling of the 
bimagnon continuum as shown in the comparison with 
the density of states (panel 6d). The correspondence of 
the Cu L 3 results with the lower ridge of bimagnon con- 
tinuum is very clear (as usual the energies in the density 
of states are in excess since the magnon-magnon interac- 
tion is neglected). This means that a particular subset 
of states is emphasized at each edge by matrix elements 
and not that some states are rigorously excluded due to a 
selection rule. In fact the calculations at the oxygen edge 
show a very tiny but finite intensity in correspondence to 
the lower ridge of the density of states. The above anal- 
ysis cannot be extended to much lower q values than in 
Fig. 6 due to the difficulty of disentangling not only bi- 
magnon from single magnon but also from phonons and 



multiphonons. In this connection we have the evidence 
that at the oxygen edge there is also a multiphonon con- 
tribution. This is not the topic of the present paper and 
for completeness is briefly summarized in Appendix B. 



B. Comparison between Cu K and Cu L 3 RIXS 

While the spectroscopy of magnetic excitations at Cu 
K has been investigated rather extensively in the last 
few years, 7,17 the comparison between L3 and K edges 
in terms of bimagnon spectral function remains an open 
issue. Indeed we have just shown that Cu L 3 explores 
the lower ridge of the bimagnon continuum; this type of 
sampling is done also at the Cu K edge accordingly to 
Ref. 19 so that a similarity between Cu K and Cu L3 is 
expected. But similarity does not mean necessarily iden- 
tity. In effect within the approximation used in Ref. 19 
the two cases should give the same spectral function apart 
a scaling factor, but this is only an approximate scheme 
and in the real world there might be differences. There- 
fore an experimental comparison between RIXS at the 
two edges is significant and is given here. In this connec- 
tion it is useful to have L3 and K spectra exactly at the 
same q values and to this purpose we measured also Cu 
K spectra. These agree with the literature in the cases 
already investigated and add new information. As it is 
well known, in the Cu K edge spectra the bimagnon is 
a small feature on top of the tail coming from the elas- 
tic peak as in Fig. 7b which gives a spectrum measured 
at the Cu K edge absorption peak (Fig. 7a). The Cu K 
RIXS spectra are given by the black open dots in Fig. 7c; 
these are raw data after background subtraction and not 
rescaled so that the relative heights of the spectra are 
significant. To these results the extracted L3 spectra are 
superimposed (open red dots); also in this case the rel- 
ative heights are preserved while the heights of the two 
sets (K and L3) are adjusted to superimpose the spectra 
at (0.6,0.6) and at (0.5,0.5) (units so that the corner of 
the BZ is at (1,1)). The main results is that the spectral 
shape is the same at the two edges within the experimen- 
tal errors dictated by statistics, which is on the average 
20%. 

It has been already demonstrated by Ellis et al. 17 that 
the accurate fitting of the bimagnon energy requires the 
inclusion of the magnon-magnon interaction so that there 
is no point in doing this exercise in the present case. 
Instead it is interesting to test at Cu K edge the basic 
approach we used in the O K case, with the obvious 
needed changes. 19,40 This approach at lowest order in 
J/T with only nearest neighbor interaction (i.e. without 
second and third neighbors and without ring exchange) 
gives the maximum of the bimagnon shown by the black 
ticks in Fig. 7. These are obtained with J— 135 meV 
from Ref. 19. As expected the calculation overestimates 
the energy at higher q. Moreover there is some under- 
estimation at low q but the general picture is remarkably 
good. This is a further, although indirect, support to 
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FIG. 7. (Color online), a) Cu K XAS of LCO. b) Typical 
RIXS spectrum with Cu K excitation at the energy marked 
in a): the gray area is the bimagnon while the blue line in- 
dicates the out of resonance elastic background measured 40 
eV below the main XAS peak, c) Cu L3 and Cu K bimagnon 
along (1,1) direction, respectively in red and black dots. The 
vertical offset between two consecutive spectra is proportional 
to the difference of the corresponding momenta (marked on 
the right). The black ticks mark the theoretical bimagnon 
peak positions obtained from theory as in Refs. 19 and 40 
with J=135 meV. 41 



the use of a conceptually equivalent approach at O K 
edge. In this analysis we refer to the peak position which 
is a safe indicator only marginally affected by the higher 
order contributions, clearly present in the measurements. 
In fact the end point of the measured spectra is around 1 
eV while the end point of the bimagnon with interaction 
and including the experimental resolution cannot exceed 
about 0.6 eV. 



mark case of undoped LCO and we have shown that the 
spectrum is given by the even order magnetic excitations 
and in particular is dominated by bimagnon. On this ba- 
sis we have pointed out the fingerprints of bimagnon exci- 
tation which are available to be used to identify bimagnon 
also in more complex systems by using O K RIXS. The 
most typical aspect is the almost complete absence of 
dispersion and another aspect is the roughly constant 
intensity with polarization vs. q (going toward grazing 
incidence) when incoming a polarized light is used. The 
behavior of the intensity is better seen by including cor- 
rections for self-absorption which however do not change 
the general picture. These typical results are understood 
by comparing directly with the non-interacting bimagnon 
density of states. This shows that O K RIXS is sampling 
the upper ridge of the density of states near the T point; 
this argument is also supported by a basic theoretical 
treatment of the type of Ref. 40. 

Moreover a more general picture of bimagnon spec- 
troscopy is obtained by comparing the O K RIXS with 
Cu edges results i.e. L 3 and K RIXS. Both Cu edges 
explore the lower ridge of the bimagnon density of states 
and show clear dispersion. In the L3 case we present orig- 
inal data obtained with a separation procedure of odd 
and even order excitations. This allows also a compari- 
son between L3 and K edge data taken to this purpose; 
in both cases the spectral functions have the same shape 
within the noise, a point never assessed before. 

The difference in the sampling of the DOS given by 
oxygen and copper sites is not due to a selection rule but 
to a drastic difference in the contributions to the spec- 
tra from the two ridges. The high energy one dominates 
the DOS and thus the spectrum as in O K RIXS, unless 
something special happens to reduce its weight. This is 
the case of copper excitation where the matrix element 
effect, as explained in Ref. 19 for Cu-K edge, almost sup- 
presses the upper ridge and in the meantime the total 
spectral weight tends to zero with decreasing q. This 
last effect comes from the commutation of the first order 
scattering operator with the ground state Hamiltonian at 
q=0. For symmetry reasons this happens also at Cu L 3 , 
but not at O K edge. 40 

This general view shows that each edge has specific 
features summarized for convenience in Table I. In par- 
ticular we stress that O K RIXS although with a limited 
momentum gives unique information. In fact a sampling 
of even spin flip excitations near Y point with energies 
around 450 meV is not possible with K and with L3 
RIXS. Thus O K RIXS has a specific role in the study 
of the low energy scale excitations in Cuprates. 
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O K 


Cu L 3 


Cu K 


transition 


Is -> 2p 


2P3/2 -> 3d 


Is — >• Ap 


typical energy 


530 eV (soft x-ray) 


930 eV (soft x-ray) 


8990 eV (hard x-ray) 


core-hole spin-orbit 


no 


yes 


no 


spin excitation 


AS=0 


AS=0,±1 


AS=0 


single magnon 


absent 


dominant 


absent 


physical information 
on the bimagnon 


upper ridge bimagnon 
(almost flat) 


lower ridge bimagnon 
(dispersing) 


lower ridge bimagnon 
(dispersing) 


momentum extension 


- 45% of the 1 st BZ 


- 80% of the 1 st BZ 


several BZs 


RIXS contribution with 
respect to total emission 


small 


large 


small 


data handling needed 
to obtain bimagnon 


straightforward 


delicate separation 
from single magnon 


delicate background 
subtraction 



TABLE I. Summary of the main aspects of RIXS at the three edges O K, Cu L3 and Cu K. 
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Appendix A: BIMAGNON EXTRACTION AT Cu 

L 3 

The aim of this Appendix is to present the concept on 
the bimagnon extraction from the spectra at Cu L3. The 
basis are the transformation properties of the spectral 
components with the angles and/or the incident polar- 
izations. This makes it possible to separate the odd spin 
symmetry (single magnon and higher order odd excita- 
tions - referred to as the S component) from the even spin 
symmetry (bimagnon and even higher order - referred to 
as the B component). 

Accordingly to Ref. 13 the magnetic intensity can be 
expressed as the product I = F x G where F is an 
atomic form factor containing the polarization effect and 
depending on the excitation one is looking for, while G 
depends only on the absolute value of the momentum \q\. 
This factorization is essential in the procedure presented 
here and is possible because the RIXS process has only a 
single intermediate state; this is an important difference 
with respect to the more general case presented in Ref. 
14. 



We consider the spectra taken at the same ab- 
solute value of the transferred momentum q and 
with the two incident linear polarizations a and 7r, 
i.e. the four cases a), (+(?, tt), (— #,cr), (— q,ir). 
Moreover we assume they are already corrected for 
self-absorption. The function G is the same for the 
quadruplet of spectra so that only the form factor F 
determines the transformation coefficients of the S and 
B components within any pair of spectra I\ and I 2 
arbitrarily extracted from the quadruplet. One can write 

h =5i + Bi, 
I 2 = S 2 + B 2 = 7S1 + f3B 1 

so that 1\ and I 2 can be obtained since the coefficients 
7 and j3 are known from the form factors given in Refs. 36 
and 42. In fact 7 = Fs 2 /Fs\. The B symmetry is shared 
by a variety of channels including the elastic, which has 
the form factor Fe defined in Ref. 13. The factor Fb is 
proportional to Fe so that /3 = Fei/Fei = Fb 2 /Fbi- 

The decomposition is safe provided that the couple I\ 
and I 2 is properly chosen among the quadruplet. A very 
convenient choice is to take I\ and I 2 at the same angle 
(i.e. the same q) and different polarization (1 = tt and 2 
= a). In this case the self-absorption is almost the same 
and this reduces the uncertainties. 



Appendix B: PHONONS AND OVERTONES 

Phonon spectroscopy with RIXS is not the topic of the 
present paper and the instrumentation was not optimized 
to this end. In spite of that a byproduct is an interesting 
hint on phonons presented in this Appendix. 

In the paper we have shown that the bimagnon con- 
tribution in the O K RIXS spectra is basically around 
3 J so that the instrumental resolving power used here 
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FIG. 8. (Color online). Decomposition of the LCO spectrum 
(red filled dot line): a) inelastic part (black open dot line) 
after subtraction of the elastic contribution (black solid line, 
FWHM=150 meV); b) further decomposition of the low en- 
ergy inelastic spectrum by subtracting the contribution from 
the phonon breathing mode at 85 meV (green line). The re- 
maining inelastic spectrum at higher energy losses (blue filled 
dot) shows the presence of further multiphonon modes (ar- 
row). The red solid line is a smoothing curve as a guide to 
the eye. 



is more than enough to study bimagnons. The situa- 
tion with phonons is quite different: the highest energy 
phonon (around 85 meV), 43 which is also the dominant 
phonon contribution, is just above the half width (150/2 
= 75 meV) of the experimental elastic peak, measured 
accurately with scattering from a carbon tape. With 
this resolution the phonon contribution is expected to be 
seen basically as an asymmetry of the peak around zero 
energy loss and dominated by the elastic contribution. 
This is indeed the case as shown in Fig. 8(a) based on 



the data at q=0. Once a symmetric elastic peak rep- 
resenting the experimental resolution is subtracted, the 
difference spectrum (open dots) shows an extra feature 
around 120 meV. We use a Gaussian elastic peak because 
this shape is given within the error bars by the carbon 
tape measurements. 

The peak of the difference spectrum belongs to the 
region comprising phonons and their overtones which 
are not separated with the present resolution. In effect 
the high energy phonon has little dispersion so that the 
biphonon is around 170 meV. A rough decomposition 
of the open dot spectrum (see Fig. 8(b)) is compatible 
with this argument since the subtraction of a peak at 
85 meV gives a reasonable biphonon signature pointed 
out by the red arrow. This exercise is resolution lim- 
ited and cannot be quantitative, but it has an important 
implication besides the obvious fact that better resolu- 
tion is needed; it shows that with about 2 or 3 times 
better resolving power the biphonon spectral weight is 
sufficient to make it observable at the oxygen edge in 
bidimensional undoped cuprates having a large superex- 
change. In fact in these systems the magnetic excitations 
are nearby the biphonons and it is not obvious a priori 
that the biphonon spectral weight is sufficient. In effect 
up to now the biphonon has been seen with RIXS at the 
oxygen edge only in systems having very low J, 44 where 
the situation is much more favorable. 

As matter of fact the needed resolution is al- 
ready available, 27 ' 28 so that the measurement of the 
biphonon/phonon ratio in bidimensional cuprates is a 
real possibility. This is extremely important in the light 
of the recent theoretical work showing that it is possible 
to extract directly the e-ph interaction matrix element 
from this ratio obtained from the RIXS spectra. 45 
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